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ABSTRACT 


In this dissertation, it will be addressed the distribution of CO, species, carbonate (HCO, -) 
and bicarbonate (CO ?), in the aqueous molecular surface for different solution acidity levels 
by x-ray photoelectron spectroscopy (XPS) using the microjet technique. The distribution of 
the species in the surface was compared with those in bulk probed by nuclear magnetic 
resonance (NMR) performed in the same samples. This is the first results obtained by the 
Group of Diluted Nanosystems (GND - Unicamp) together with other researchers. 


For the first time, five samples of sodium carbonate 0.5M at different pH are studied. In 
the XPS C1s spectra were observed three features belonging to three different protonated 
species of CO, molecules in the solution. The chemical shift between the carbonate and 
bicarbonate molecules were enough to produce an unambiguous identification of the 
Cis areas and thus their amount on the molecular surface. The analysis of XPS spectra 
showed the variation on the concentration of the molecules as a function of solution acidity. 
Comparing with bulk the relative concentrations presented clear differences, especially for 
higher pH. The data were fitted using an enrichment model proposed in the literature. From 
this fitting, it was verified that doubly charged species are less prone to stay on the surface 
compared to the singly charged in agreement with previous theoretical predictions for ionic 
salts. It was proposed a new procedure to obtain the Gibbs free energy for two of the 
protonated species. This new model is a variation of the well known Langmuir model and 
results obtained here are consistent with the enrichment model and may be of general use. 


The distribution of the same species in bulk was studied using NMR. The samples used 
are the same as the XPS experiment. Five measurements of carbon 13 were performed at 
room temperature. Based on the chemical shifts, it was possible to obtain the concentration 
of the species as a function of the bulk acidity. The Hendersen-Hasselbalch equation which 
describes the bulk acid-base reactions was chosen for the data analysis. 


The observed differences between bulk and surface are discussed including a possible 
higher uptake of CO, by aqueous solutions as previously thought. This finding will obviously 
have consequences for vital systems such as oceans or atmospheric aerosols. This new 
revealed scenario may increase awareness of the urgent measures needed to avoid some 
of the deleterious impacts of acid rains or acidification of oceans sooner than current models 
predict which are based on concentrations found in the solution's bulk. 


KEYWORDS: PES, photoelectron spectroscopy, XPS, liquid surface, microjet, carbonate, 
bicarbonate, carbon dioxide, deprotonation, pH, pK,, liquid bulk, Langmuir. 


RESUMO 


Nesta dissertação será abordada a distribuição de espécies de CO,, carbonato e 
bicarbonato, na superfície molecular aquosa para diferentes níveis de acidez do meio 
através da espectroscopia de fotoelétron (XPS) usando a técnica de microjato. Para a 
comparação com a distribuição destas espécies no corpo da solução um experimento de 
ressonância magnética nuclear (RNM) foi realizado nas mesmas amostras do experimento 
anterior. Estes são os primeiros resultados obtidos pelo Grupo de Nanossistemas Diluídos 
da Universidade Estadual de Campinas (GND - Unicamp) em conjunto com outros 
pesquisadores. 


Pela primeira vez foi realizado um estudo das moléculas de bicarbonato e carbonato, 
provenientes da reação de deprotonação do ácido carbônico na superfície aquosa. Foram 
preparadas 5 amostras de carbonato de sódio a 0.5M em diferentes pHs que foram 
utilizadas para realizar o experimento de XPS. Nestes 5 espectros de C1s foi observado a 
presença das três moléculas na solução, sendo o deslocamento químico entre carbonato e 
bicarbonato suficiente para fazer uma identificação de cada pico no espectro. A análise dos 
espectros de XPS mostrou a variação da concentração das moléculas em função da acidez 
da solução. Foram observados diferenças significativas principalmente para pHs mais altos 
em relação aquelas encontradas para o interior do líquido. Os dados foram ajustados com 
o modelo de enriquecimento da superfície, proposto em um trabalho na literatura. A partir 
deste ajuste, foi constatado que a espécie duplamente deprotonada é menos propensa a 
se concentrar na superfície do que a espécie com deprotonação simples. Ambas tem sua 
concentração reduzida na superfície em relação ao interior da solução. Os dados também 
foram ajustados com um novo modelo proposto baseado no modelo de Langmuir. A partir 
deste ajuste é possível obter a energia livre de Gibbs na superfície que se comparou 
positivamente com o modelo de enriquecimento. Esta nova forma de tratamento dos dados 
pode ser aplicada a outras moléculas sendo portanto de interesse mais geral. 


Também foi realizado um estudo sobre a distribuição das mesmas espécies no corpo do 
líquido usando RMN, que uma técnica experimental sensível ao interior da solução. As 
amostras utilizadas são as mesmas do experimento de XPS sendo que cinco medidas 
de carbono 13 foram feitas a temperatura ambiente. Com base no deslocamento qumíco 
dos picos é possível obter a concentração das espécies em função da acidez do corpo 
da solução. A equação de Hendersen Hasselbalch que descreve o comportamento das 
reações de ácido e base para o corpo da solução foi escolhida para fazer o ajuste destes 
dados. 


Por fim, analisando as diferenças encontradas entre os resultados experimentais e 
seus respectivos ajustes dos dois experimentos foi constatado que as duas moléculas 
apresentam concentrações e dinâmicas diferentes na superfície e no interior. 


Com bases nos resultados obtidos foi discutido na tese que é possível prever que a 
adsorção de CO, em soluções aquosas tende a ser maior do que aquela prevista por 
modelos baseados na concentração do interior da solução. Este fato trás consequências 
em relação à urgência de se acelerar ações para mitigar o aumento da presença de CO, 
na atmosfera tendo em vista reduzir o aparecimento da chuva ácida e acidificação dos 
oceanos. 


PALAVRAS-CHAVE: PES, espectroscopia de fotoelétron, XPS, superfície líquida, microjato, 
carbonato, bicarbonato, dióxido de carbono, deprotonação, pH, pK,, bulk líquido, Langmuir. 


CHAPTER 1 - INTRODUCTION 


It is widely accepted that the Industrial Revolution brought to humanity a great 
scientific and technological progress. Cheap and abundant fuel, such as charcoal, powered 
the steam engines and trains. The model for energy production at that time requires the 
burning of charcoal which, among other pollutants, releases carbon dioxide (CO,) in the 
atmosphere. This process has intensified enormously after the discovery of oil field over 
the years. The main source of greenhouse effect is the resulting CO, released to the 
atmosphere [7]. The concentration of atmospheric carbon dioxide is almost 40% higher 
than those found in the pre-industrial period. A significant part of this CO, added in the 
atmosphere is also absorbed by the ocean and land biosphere [8, 9]. Indeed the oceans 
are responsible for the uptake of as much as 27% of the anthropogenic CO, emissions 
[10] if a time span of a decade is taken into account. The importance of the ocean for CO, 
absorption on a millennium time scale is substantially increase since the CO, equilibrates 
largely between atmosphere and ocean. In fact, typically about 20 per cent of the added CO, 
to the atmosphere stays there and as much as 80 per cent end up in the ocean [11]. The 
North Atlantic Ocean presents the higher concentrations of CO, [12]. 


The dissolution of CO, in the ocean modify the seawater chemistry by the 
increase carbon acid formation, and as a consequence, the pH of the water decrease. 
This water becomes more corrosive to organisms that produce shells, see figure 1, and 
others calcareous structures. The low pH also affects the reproduction, physiology and 
geographical distribution of marine life [13, 14, 15]. lf the emission of CO, is maintained, in a 
few decades the sea water in tropical regions will not sustain the development of coral reef 
ecosystems any more[16, 17, 18]. The decrease in the ocean pH will lead to biodiversity 
loss and economic impacts. The pH of the ocean surface is 0.1 pH units smaller since the 
beginning of the Industrial Revolution, which represents a 30% change in the ocean acidity. 
Itis expected a drop of 0.3 or 0.5 pH units by the end of the century [19]. 


The absorption of CO, by the ocean reduces the impact of the greenhouse gas to the 
climate changes. However, these mitigation actions on the climate do not have any effect on 
the ocean acidification (OA). The OA is a problem as severe as the greenhouse effect which 
needs to be more deeply studied. To understand the conditions in which the CO, absorption 
occurs it is necessary to access what is present in the ocean surface on a molecular scale. 
More precisely it is desirable to determine how carbon dioxide and its protonated species 
are distributed on the molecular surface of aqueous model-solutions as a function of the 
pH. So far we assumed that these distributions are the same in bulk and at the molecular 
surface. This assumption needs to be tested experimentally since the uptake of CO, is a 


surface phenomenon. 
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Figure 1.1: This is the pteropod, or sea butterfly. Figure shows what happen to a pteropods 
shell when de sea water is in the conditions predicted for the year 2100. The acidity slower 
dissolves de shells over 45 days. Photo by David Littschwager [1]. 


Among several experimental methods commontly utilized in the study of the dynamical 
and structural properties of aqueous solutions, photoelectron spectroscopy (PES) is a 
powerful technique due to its element specificity as well as surface sensitivity. 


During the last decades, PES has its application restricted to solids and gases due 
to its high vacuum demands on the experimental chamber. Many liquids have a high vapor 
pressure at room temperature, in particular, water vapor pressure at room temperature is 
as high as 22.3 mbar [20]. At these conditions, the photoelectrons have a very short mean 
free path turning PES studies impractical. The most prominent advance to allow PES in 
aqueous solutions was the development of micro-jet technique by Faubel [21]. In this thesis, 
the surface sensitivity of PES applied to liquids will be used to study the distribution of CO, 
protonated species at the aqueous surface. This distribution willbe compared to that obtained 
from nuclear magnetic resonance spectroscopy (NMR) which is bulk representative. The 
samples are the same in both techniques. 


Chapter 2 presents a brief discussion about water and carbon dioxide electronic 
structure as well as some basic concepts needed in this work such as protonation and 
de-protonation processes in aqueous solutions. Chapter 3 presents a review on the PES 
technique and a few important concepts related to the surface sensitivity and the element 
specificity of this method. Chapter 4 presents a brief history of the development of micro-jet, 
the type of photon source used in this work, the conditions of the liquid surface atthe moment 
the measure is carried out, the curve fitting process employed and the model used to fit the 
PES data and finally there is a brief review of the NMR concepts, the experiment and the 
fitting procedure adopted. Chapter 5 presents the results obtained with both experimental 
techniques and the discussion of the results. The conclusion and some perspectives are 
shown in chapter 6. 
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CHAPTER 2- BASIC CONCEPTS 


This chapter deals with basic chemical-physics concepts of liquid water. Initially, a 
description of the water molecule and the liquid phase are given, then a description of 
CO, molecule, next some important chemical concepts and a brief discussion about the 
deprotonation processes of carbon dioxide in aqueous solution. Finally, we introduce some 
concepts required to describe the liquid interface. 


WATER 


Matter can assume the most different states of matter such as solid, liquid, gas 
and plasma. Some others states observed in extreme situations such as Bose-Einstein 
condensates, neutrondegenerate matter and quark-gluon plasma. For this dissertation it 
will be considerate only the first four states of matter which are the most observable in the 
every day life. In solids, the interaction between particles is strong, not allowing too much 
mobility unless in the sites of the crystal lattice. Gases particles almost do not interact (of the 
kg T order) and move independently. Liquids are hybrids of solids and gases. The interaction 
energy between particles is comparable with thermal energy like gases but with density 
closer to solids. 


Water is liquid at room temperature, odorless and known as the universal solvent. It 
is a planar molecule with C, symmetry, composed of two hydrogens attached to an oxygen 
atom by covalent bonds. Those bonds form an angle of 104.5Ainthe gas phase and a length 
of 0.957 85 A [22]. In liquid phase those values are slightly higher as shown by ab initio 
calculus and diffraction studies [23, 24, 25]. Due to the oxygen electronegativity, there is a 
charge shift resulting in the presence of a permanent dipole moment in the water molecule. 
The molecular symmetry defines its orbitals form. The water electronic configuration in the 
ground state is 


(la) (2a P(1b) (3a, (1b,)? (2.1.1) 


The orbital 1a, does not participate in the chemical bonds and have an atomic 
characteristic; it can be referred to as the O1s water orbital. The orbitals 2a, and 3a, are 
partially binding and nonbonding. The orbital 1b, is stricily binding and the valence orbital 
1b, has characteristic of the O2p, orbital, known as lone pair and is nonbonding. 


Despite the apparent simplicity of water molecules, liquid water gets more 
complicated due to the formation of hydrogen bonds. The formation of these bons occurs 
when a partially positive hydrogen from one molecule is positioned towards the lone pair 
electrons belonging to the oxygen from another molecule. Each water molecule can form 
four hydrogen bonds, two as donor and two as receptor. Due to the linear geometry of the 
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interactions and molecular geometry the final structure that only one molecule does four 
bonds is a tetrahedral which produces the crystalline form of ice, as shown in figure 2.2. 
In liquid water, the tetrahedral structure is locally found, and the increasing temperature 
reduce the formation of this structure. 


(0) H20 H-H 


Hz às 


Ois laí 


Figure 2.1: Diagram of energy levels of occupied molecular orbitals of water gas phase. 
The linear combination of hydrogen and oxygen atomic orbitals formed the molecular water 
orbitals. 
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Figure 2.2: (a) Model of tetrahedral structure of hydrogen bonds between water molecules. 


(b) Real-space visualization of the formation of hydrogen bonding in 8-hydroxyquiline (8-hg) 
molecular assemblies on a Cu(111) substrate using noncontact atomic force microscopy (NC- 
AFM) [2]. 
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The hydrogen bonds also occur when the hydrogen atom makes covalent bounds 
to nitrogen or fluorine, which are atoms highly electronegative. The enthalpy of hydrogen 
bonds is between the covalent bonds and the Van der Waals [26]. This bond occurs in 
organic molecules, which plays a vital role in secondary and a tertiary structure of DNA 
and proteins molecules [27], and inorganic molecules, which is responsible for the highest 
boiling point of water compared to other liquids [28]. The size of hydrogen bonds varies 
depending on the number of bonds a molecule makes; it is typically given by 2.8 À, which 
is shorter than a Van der Waals bond length of 3.65A [26]. 


CARBON DIOXIDE 


The carbon dioxide (CO,) is a linear molecule with symmetry D., very 
thermodynamically stable and nonpolar. Itis composed by two atoms of oxygen attached 
to one carbon atom by covalent bonds. It is an odorless and colorless gas present in the 
atmosphere, anditis an important component of plants and animals live cycle. The process of 
cellular respiration and by the combustion of fuels in particular of fossil origin produces CO. 
Carbon dioxide is solid at -78.5'C and sublimates 1 atm being widely used in the industrial 
process that requires largescale refrigeration, in this state, it is called dry ice because it 
sublimates at room temperature. Its density is higher than the air, so it accumulates in 
environments without ventilation. The solution of CO, in water will be better described in 
section 2.4. The double bonds C = O of the molecule are equivalent and with a 116.3 pm 
length [29]. Carbon dioxide ground state configuration is: 


(19 (10,)(20920, MIM) (2.2.1) 


As said before, the linear combination of atomic orbitals form the molecular orbitals. 
Figure 2.3 presentes an estimation of the molecular orbital energy levels diagram of the 
CO, molecule. The first ionization potential of carbon dioxide is 13.778 + 0.002 eV [30] so 
this molecule is a rather weaker electron donor. The Lowest Unoccupied Molecular Orbital 
(LUMO) is 3.8 eV which is a rather low value translating in a high electron affinity which is 
associated with the central carbon atom. Thus CO, in a reaction with transition metals is 
considered a poor electron donor but a good electron acceptor [31]. 
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Figure 2.3: Molecular orbital energy levels diagram of occupied energy levels in the gas 
phase of CO,. 


For the carbon atom to perform four bonds with two oxygen atoms and to keep the 
linearity of the molecule the carbon atomic orbitals need to be hybridized, that is, the 
mixed orbitals form a new hybrid orbital. This mixture is a linear combination of the original 
orbitals. The ground state of carbon is 15252p2p,, see figure 2.4. Energy is necessary 
to promote one electron from 2s orbital to 2p orbital. The only orbital available to form the 
hybrid orbital is p, forming the sp orbital, once p, and p, are fullly occupied. Thus the carbon 
atom has four electrons paired, and now it can do the four bonds required to the formation 
of the CO, molecule, the two bonds made by the pure orbitals are mi bond and the other two 
bonds made by the hybrid orbitals are o bonds. The hybridization of the carbon atom is 
energetically favorable because it minimizes the energy of the molecule. 
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Figure 2.4: Hybridization of carbon atom to form the CO, molecule. First is the ground 
state and the electron distribution in the orbitals without hybridization. Below is the electron 
distribution in the orbitals with the hybridization. 


CHEMICAL CONCEPTS 


In this section introduces some concepts used by chemists as concentration, the 
solubility of a solute and chemical equilibrium that will be largely used in this thesis. 


Molar concentration or molarity (M) is the ratio between the amount of solute by the 
volume of the solution [mol/dm*][32]. This concept is easily related to results in photoelectron 
spectroscopy. On the same line there is the molar fraction that is a dimensionless number 
given by the ratio of solute moles by the sum of the solute and solvent, see equation 2.3.1 
[33]. The concentration does not depend on temperature. 


Ni 


X=— i=1,2 
ny + nao E (2.3.1) 


the molar fraction can be calculated on the solute with i = 1 or on the solvent with i = 2. 


Solubility is the maximum amount of a substance dissolved in liquid [32]. The polarity 
determines if a substance will be more or less dissolved in a solvent. Polar substance 
dissolves better in polar solvents and apolar substance dissolves better in apolar solvents. 
The solubility depends on temperature and others dissolved substances. The saturation 
happens when solute has the maximum concentration, and the solution can not dissolve 
any more solute. 


Gases also are dissolved in liquids. Gases molecules are in constant movement, 
some of those molecules with enough velocity dissolve by penetrating the liquid surface. 
This mechanism depends on the gas pressure (the higher the pressure the greater the 
strength of the collisions), the liquid temperature (higher temperature induces more 
movement of molecules which increase the chances that gas molecules escape the liquid) 
and the fluid surface agitation which may facilitate the gas exchange. The gas dissolution in 
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liquid is governed by Henry's law [32], equation 2.3.2. In this law, the gas solubility depends 
on the gas partial pressure exerted on the liquid. The constant of proportionality k,, Henry's 
constant, varies with the gas type and temperature. 


Co =k 


sol H , Pas 


(2.3.2) 

The combination of reagents to form a new substance with different properties is 
called chemical reaction [32]. Some reactions are reversible, in which case the reagents 
combine to create products, and the products combine to regenerates the reagents. It is 
a dynamic process; the reaction happens in both ways. Each reaction has its velocity, and 
it depends on the concentration of reagents. The chemical equilibrium is achieved when 
forward and reverse reactions have the same speeds. Considering the hypothetical reaction 


A+B=C+D (2.3.3) 
thus the equilibrium constant is given by 


p= AIIBI 
CID] (2.3.4) 


where [A], [B], [C], [D] are the molar concentrations. This constant is temperature 
dependent, but do not depend on initial reaction conditions. Greater the equilibrium constant 
value, higher is the formation of products and higher is the reaction yield. 


In Brônsted-Lowry theory [34], acids are defined as proton donors and bases are 
defined as proton acceptors which can ionize when dissolved in water. The ionic conjugated 
formed in an aqueous environment is stable in contrast with gaseous phase. Considering 
a generic dissociation reaction of an acid in water forming H,0*, A- and its equilibrium 
constant 

[H50*][A7] 


= Hot — É = 
HA+H0S HO! +A" and Ki HA] (2.3.5) 


The dissociation of HA originates equal amountis of H,O* and A”, thus they have the 
same concentration, [H,0*] = [A]. Applying the logarithm in equilibrium constant equation 
and using that log[H,0*] = pH and logK, = pK, the famous Hendersen-Hasselbalch equation 
[35] is obtained. 


HA 
pH = pK, + log (55) (2.3.6) 
Another way to study the relation between the acid and base concentration as a 
function of the sample's pH is through the acid fraction, which is the ratio between the acid 
concentration by the total concentration, as 


[HA] 


E= TAF 


(2.3.7) 
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Replacing equation 2.3.6 in the equation above it is found a relation between the acid 
fraction as a function of the pH on the sample. 


1 


= 10PH-pKa SA (2.3.9) 


DEPROTONATION 


Deprotonation is a chemical term that refers to the removal of one proton from an 
acid in an acid-base reaction forming a conjugated base. Protonation is the process of a 
proton added to base forming the conjugated acid. Species that can both accept and donate 
a proton are called amphoteric, as the water molecule that receives a proton forming the 
hydronium ion H,O* or lose a proton forming the hydroxide ion OH-. By the Henderson- 
Hasselbalch equation 2.3.6 when the acid concentration and your conjugated base are 
equal then we have pK = pK,. When gaseous carbon dioxide is dissolved in water a rapid 


equilibrium occurs CO, = CO, q? then this dissolved carbon dioxide and water react to 


2(9) 
form the carbonic acid (H,CO,) in chemical equilibrium. 


CO, ay + HO = HCO, (2.4.1) 

which equilibrium constant is Ka =1.7 10% in pure water [36]. The small value of that 
constant indicates that the majority of the carbon dioxide is not converted to carbonic acid. 
Carbon dioxide is a polyprotic acid, in other words, it can donate more than one hydrogen 
ion forming bicarbonate (HCO) or carbonate (COZ). The relative concentration between 


those species depends of solution pH as is shown in figure 2.5. 


Fractions 


Figure 2.5: Variation of species concentration H,CO,, HCO; and COZ as a function of pH 
liquid bulk. 
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The deprotonation reaction of carbonic acid is presented at figure 2.6 below. Each 
reaction has its own dissociation constant, for the first reaction k, = 4.45 x 10-/mol.l and 
for the second reaction k, = 4.84 10-"'mol.l” [29]. The k, is larger than k, because for the 
second deprotonation the proton is separated from a center that is more negative than in the 
first deprotonation, so this additional work makes the second deprotonation less favorable. 
Rather than writing the value of the dissociation constant, it is prefer to use the definition 
pK, = logk,. The pK,, value is usually quoted as 6.37, this evaluation assumes that all the 
acid present in solution is carbonic acid but, as said before, there is a large proportion of 
dissolved CO,. Taking this in account the correci value is pK,, = 3.6. The pK,, value is 10.3. 


carbonic acid bicarbonate carbonate 
O « 
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Figure 2.6: Deprotonation reaction of carbonic acid 


Van't Hoff equation 


Those pK, values are dependent on temperature, which dependency is given by the 
Van? Hoff equation. But before deep in the consequences of this equation, it is important to 
do a brief review on some concepts that concern chemical equilibrium. This review can be 
found in any chemistry book, such as in reference [32] that is more detailed, but here it is 
generally shown the concepts to help the reader that is not acquainted with the subject to 
understand the experimental results presented at the end of this dissertation. 


The minimal Gibbs energy gives the chemical equilibrium of a reaction. Considering 
a common equilibrium reaction as A = B, after some time an amount of A will turn into B, 
so the amount of each one will turn to dn, = -dã and dn, = +d&, where & is the extended 
reaction. The slope of the Gibbs energy as a function of extended reaction is the reaction 


OG 
A,G — Es 
( os ). RA] 


The above equation is valid for constant pressure p and temperature T. The total 


Gibbs energy. 


derivative of Gibbs energy is given by 
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dG = u,dn, + u,dn, = (4, — uo)dE (2.4.3) 


where uy, is the chemical potential for reagents and products. From the equation 
above one can identify 


pi (2.4.4) 

The equation 2.4.4 allows a different interpretation of reaction Gibbs energy, which 
is the difference between the chemical potentials. As the reaction proceeds the chemical 
potentials change as well as the reaction Gibbs energy, the reaction tends in directions 
that decrease G. The chemical equilibrium is achieved when the AG = 0, or when u, =H,. 
In other words, the chemical equilibrium is reached when the chemical potentials are in 
balance. 


IfA and B are perfects gases the chemical potential can be written as y = uº + RTInp, 
so the reaction Gibbs energy is 


AG = (4º, + RTInp,) — (u,º + RTInp,) = A,Gº + RTInQ (2.4.5) 


where Q = p,/p, is the reaction quotient of partial pressures p. Despite the 
approximation with perfect gases, the final result can be obtained from an general reaction. 
In the chemical equilibrium AG = 0 so 


AGº = -RTInK (2.4.6) 


the K constant is the ratio between the partial pressures when the reaction is at 
equilibrium. Equation 2.4.6 is an important thermodynamic relation between the reaction 
Gibbs energy and the dissociation constant. 


Changes in temperature, pressure and the concentration of reagents and products 
may have some effect on the chemical equilibrium and in the dissociation constant. The 
chemical equilibrium was defined as a ratio between the partial pressure of the reagents and 
products, so a pressure variation may tend the reaction to a production of more reagents 
or products, but the ratio between the partial pressure will stay the same, therefore K is 
independent of pressure variation, the same way with A or B concentration, which is a 
consequence of the Le Chatelier's principle. 
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Temperature has a different effect on dissociation constant. The addition or the 
removal of energy disturbs the reaction, and from the Le Chateliers principle, the reaction 
is opposed to this disturbance of equilibrium. In the case of an exothermic reaction, the 
increase of temperature will favor the reagents production, because the excess of energy is 
absorbed as heat to oppose the disturb occurs by the rising temperature. For endothermic 
reactions, the same increase in temperature cause a favor in products production. To 
understand this effect 2.4.6 is derived to get 


dnkK — 1MAG/T) 


dT Rar pe 
and using the Gibbs-Helmholtz equation 
0AG/T AH 
oT “q (2.4.8) 
p 
the Van't Hoff equation is given 
dink A, Hº 
d” — RT? a 


From this equation it can be seen that 


* For exothermic reaction, A, 1º <0, the &É < 0, which means that InK decrea- 
ses when the temperature increase. 


* | For endothermic reaction, A,Hº > 0, the &“* > 0, which means that InK increa- 
ses when the temperature decrease. 


The constant dissociation k, at temperature T, in terms of and initial conditions k, and 
T value is obtained by integratind equation 2.4.9. 


1/To 


1 
Inks == Ink = "“R i A,Hºd(1/T) (2.4.10) 
á 1. N 


assuming that A Hº is almost constant over the temperature range of interest 
1/T» 


1 
Ink, — Ink = E A Hed) (2.4.11) 
UT 


In this dissertation, the studied reaction is the deprotonation reaction of bicarbonate 
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into carbonate in bulk and liquid surface, so XPS experiment was performed to observe the 
distribution of both species in the liquid surface and a NMR experiment that is bulk sensitive. 
The experiments were performed at different temperatures so it might have some influence 
on the pK, value. The van't Hoff equation will be used to predict the pK, value and then 
it will be compared to that obtained from the data fitting. It is first necessary to determine 
if the deprotonation reaction is exothermic or endothermic to predict the behavior of the 
dissociation constant at different temperatures. A study shows that the temperature increase 
gives an increase in the dissociation constant [37], as previously mentioned, reactions that 
present this behavior are given by endothermic reactions. In the same study, the enthalpy of 
the dissociation reaction was calculated using vant Hoff equation, the value found was A,H 
= 13.7 + 0.5kJ.mol”", closer to the value from others studies [38]. 


LIQUID INTERFACE 


Before some concepts about the vapor-liquid interface are introduced it is important 
to delimit the meaning of surface and bulk, that will be largely used in this thesis. Considering 
a liquid water jet with radius R in vacuum, the density p(r) is described as follows: in r = O 
the density is equal to liquid water and in r > R the density is equal to the water vapor. In 
this case, under experimental conditions, the vapor pressure varies proportionally to 1/r in 
this region. The region at r = R is the vapor-liquid interface, that is the region called liquid 
surface. Bulk is the region below the liquid surface, but there is no consensus about how 
deep the atoms begin to be part of the bulk, it depends on the liquid structure. 


As said before liquid water is a very complex system, the description of bulk structure 
and interactions with ions and other solutes are complicated, however, the surface-liquid 
interface is even more challenging. The forces acting into the species residing closer to the 
interface are spatially asymmetric, unlike what happens in bulk in which the interactions of 
one particle spherically symmetrical are identical in all directions, so the resulting net force 
is null. The variation in free energy as a function of position relative to the surface is called 
mean force surface potential. As a consequence of this energy variation, some places will 
be more energetically favorable than others generating a density variation of solutes near 
the interface. The species that have a high propensity to reside in the surface are called 
surface-active. They accumulate in interface decreasing the surface tension compared with 
pure water. On the other hand, some others species have the propensity to move to the 
bulk. This movement of molecules is called adsorption. An obvious consequence of this 
movement is an altered spatial distribution of solutes at the interface. 
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Figure 2.7: Liquid jet in a vacuum with radius R, which is the interface that separates the 
liquid region from the vapor. 


Considering the liquid interface as an infinity planar surface and dimensions x and 
y are perpendicular, the relation between the mean force surface potential w(z) and the 
density profile n(z) of the |" solute with bulk concentration n, at an Z distance from the 
surface is given by 


ni(z) = 4 - exp (40) 


w, governs the surface propensity, but it depends on the nature of particle and 


(2.5.1) 


interface and it varies with solute concentration. 


There are two concepts of surface tension. The mechanical one says that surface 
tension is the force applied to a surface plane per unit length. The thermodynamic one is 
defined as the excess of free energy due to the surface, which entails ffom the attraction 
and repulsion that molecules in the surface feel. Combining these two concepts, the surface 
tension o is expressed as the work necessary to change the surface area o = dW/dA. The 
surface tension changes if the concentration of species is not uniform throughout the liquid, 
there may be an accumulation or lack at the interface. This variation can be related to 
chemical potential by Gibbs isotherm [32] 


do 
oo e a (2.5.2) 
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Where T, is the amount of solute per unit surface area. This equation relates the 
surface tension between a solution with the solute surface excess T;, and the chemical 
potential of the i" component y, [39] which implies that the surface tension decreases when 
the surface excess concentration is positive, and it increases when the surface excess 
concentration is negative. 


So far the description of the surface is very rough, a better description of the 
microscopic surface structure is needed and then one can find a density profile expression 
of ions that is closer to the interface. 


The electrostatic model published by Onsager and Samaras [40] is based on the 
classic image charge problem. The ions at interface are charged points, and the interface 
itself is treated as an interface that separates two dielectric media. The ion interacts with your 
image charge at the other side of the interface, and the force generated by this interaction is 
the only one that contributes to the surface potential. Following this argument, they derived 
the analytic expression below to the density profile 


qexp(—2Rz) 
l6recokBTz 


RR, 


(2.5.3) 

where q is the ions charge, k in the exponential damping factor and it comes from 
the ionic screening in solution, which comes from the Debye-Huckel theory and the factor e 
is the water dielectric constant. 


0 5 10 20 25 30 


15 
Figure 2.8: Superficial density profile for two concentrations of one Debey-Huckel 


electrolyte according to the Onsager-Samaras prevision to 300K. Note that on the surface 
(at z = 0) there are no ions. 
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Figure 2.8 is the representation of equation 2.5.3. In the interface, the density 
decreases according to the Gibbs adsorption equation 2.5.2. This model describes very 
well the general behavior of surface tension in aqueous solutions, but for very concentrated 
solutions the model does not fit correctly. In this dissertation, the studied solutions present 
low concentration so the model can be applied to this case. 
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CHAPTER 3 - PHOTOELECTRON SPECTROSCOPY AND LIQUID 
SURFACE 


In this chapter, the background of the spectroscopic techniques used in this thesis 
will be discussed. The results presented here are based on electron spectroscopy involving 
core-level ionization (C1s) while normalizations and energy calibration relied upon valence 
electron ionization. All performed on liquid samples. The photon source used in all cases is 
synchrotron radiation. The studied solutions were further scrutinized by nuclear magnetic 
resonance (NMR). NMR is only bulk sensitive therefore this technique provided a way to 
monitor the situation at that part of the solution. Firstly a review of the interaction between 
light and matter is presented mainly focused on the photoelectric effect which is the 
phenomena behind X-ray Photoelectron Spectroscopy (XPS). Particular attention is given 
to accurately characterize the surface sensitivity of XPS a well as its chemical element 
specificity. The shape of the photoemission lines is also discussed in this section. 


ABSORPTION OF LIGHT BY MATTER 


XPS is based on the interaction between light and matter and constitutes an important 
technique used by scientists to study the second. One way that light interacts with matter is 
by absorption. When an atom absorbs a photon, an electron may be promoted to an excited 
discrete state or the continuum. If the atom is excited after a period it will decay back into 
the ground state. This relaxation may occur by different paths, see figure 3.1. From the 
knowledge of the different ways of relaxation important information about the system is 
obtained. The process of absorption depends on the photon energy and binding energy of 
the electrons involved in the excitation. In case the incident photon energy is larger than the 
electron binding energy and ionization will occur which is the phenomena describing XPS. 
The ejected electron will have a kinetic energy K, given by: 


K=hv-B (3.1.1) 


where hv is the photon energy and B, the electron binding energy. This is the well 
known Einstein's photoelectric effect equation [41]. 


In general low photon energy, i. e. 10-40 eV, is used to maximize valence electron 
ionization. Valence electrons are delocalized and participate in chemical bonds. In case 
the photoionization takes place in an electron belonging to the highest occupied molecular 
orbital (HOMO) further relaxation of the final ionic state is not possible. A different history 
occurs if an inner valence electron is ionized. Further decay processes such as ionization 
happen. 
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Figure 3.1: Various ionization processes from its ground state with charge q by x-ray photon 
absorption. The ejection of an electron increases the charge of the system by +1. 


For higher photon energies inner shell electrons may also be ionized. This process, 
called core-ionization, is often used to perform “chemical analysis”. Indeed by determining 
the kinetic energy of the ejected photoelectron, its binding energy can be easily obtained 
provided that the incident photon energy is known. The core-electron binding energy is 
significantly different ffom neighboring atoms in the periodic table thus easily allowing 
element specific identification. Core-electrons are also localized in a particular atom but 
despite this localization and consequent negligible participation in the chemical bonds 
their binding energy varies according to the chemical environment. This energy variation 
is called chemical shift. A more detailed explanation of the chemical shift is presented in 
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section 3.3, but it is convenient to introduce an example of this effect. Figure 3.2 shows the 
Nis and C1s photoelectron spectra of aqueous glycine at pH 1, 7 and 13, where each pH 
present a different deprotonation state of the molecule as explained in section 2.4. Glycine 
has two carbons in its composition identified in the Cts spectrum as the two peaks. The 
deprotonation of one of the carbons causes a shift in the binding energies that is observed 
in the spectra. 
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Figure 3.2: (a) Aqueous glycine N1s and C1s photoelectron spectra at pH 1,7 and 13!. 
Reproduced from: On the origins of core-electron chemical shifts of small biomolecules in 
aqueous solution: insights from photoemission and ab initio calculations of glycineag. Figure 

reproduced from reference [3]. (b) Deprotonation process of glycine for different pH. 


The core-hole lifetime is usually a few femtoseconds according to the uncertainty 
principle, see section 3.4. After this time usually a valence electron fills the hole and the 
excess energy is released through the emission of a second electron, this process is called 
normal Auger decay [42, 43], if however a photon is emitted in this second relaxation process 
we call it normal x-ray emission or communaly named fluorescence. 


Another essential relaxation process in the presence of a formed core-hole is when a 
valence electron belonging to a neighboring atom or molecule participates in the relaxation 
process which is known as Interatomic Coulombic Decay ICD [44]. This process was 
observed recently in water [45, 46]. 


1.Zwitterion is a neutral chemical compound, but it has opposite charges in different atoms. Usually the term is used 
when the charges are at non-adjacent atoms. 
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SURFACE SENSITIVITY 


When the X-ray photon is eventually absorbed, the created photoelectron may 
undergo elastic and inelastic collisions in its way to the liquid surface. As discussed before 
the kinetic energy of emitted photoelectron is essential information that is altered by inelastic 
scattering. In many cases elastic and inelastic scattering will prevent the photoelectron from 
reaching the electron analyzer. In opposition to this situation, photoelectrons ejected close 
to the surface will have a low probability to be scattered. 


X-rays can penetrate relatively long distances inside liquid samples without inelastic 
collisions, whereas electron needs a higher energy to be able to travel the same distance. 
Figure 3.3 shows that a photon of 520 eVs can penetrate a few micrometers of a sample, 
while electrons need a 400 keVs to travel 0.5 um before being scattered inelastically. 
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Figure 3.3: Penetration depth of x-rays and electrons in water and protein. Electrons with KeV 
energy penetrate only 1 micron thickness because of multiple scattering. Reproduced from 
reference [4]. 


Two concepts are used to express where the created photoelectrons come from, 
an idea often called surface sensitivity. The effective attenuation length (EAL) means the 
shortest distance between two points that the intensity of the signal is reduced by 1/e [47] (e 
is the Euler number). The inelastic mean free path (IMFP) is defined as the mean distance 
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photoelectrons travel through a material before an inelastic collision takes place. If elastic 
scattering is efficient in the medium, the IMFP may be higher than EAL. But in the case 
of weak inelastic scattering, IMPF and EAL are almost identical as shown in figure 3.4a. 
The experimental distinction between both concepts in total signal attenuation is virtually 
impossible. Therefore the knowledge of EAL is sufficient to determine the degree of surface 
sensitivity. 

The photoelectron intensity varies as a function of the electron kinetic energy for a 
single system component and it is used to extract relative EAL values. Photoemission signal 
can be expressed by being proportional to the atomic density exponentially attenuated with 
the effective attenuation length. In the case of a plane interface, the signal expression is 
given by 3.2.1. 


I= É 2exp (-—2— 
DT) P( Jezp ( a] (3.2.1) 


Where p(z) is the density profile and a is a constant that depends on the differential 
photoionization cross section, the total photon flux, the experimental alignment and the 
transmission function of the electron analyzer [5]. For different materials, EAL has a very 
similar depen dence on electron kinetic energy, this dependence is known as the universal 
curve 3.4b, where the escape depth is minimal (5 - 10 “A) when the photoelectron kinetic 
energy is between 50 - 100 eVs, resulting in a high surface contribution to the signal. For 
higher energies the EAL increases and photoelectrons from bulk begin to contribute to the 
signal. 
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Figure 3.4: (a) Effective attenuation length and inelastic mean free path in a medium. I(A) is 
the intensity in point A and I(B) is the intensity in point B. (b) Effective attenuation length as a 
function of electron kinetic energy for water, aluminium, polymethylmethacrylate (PMMA) and 

the universal curve. Logarithmic energy scale. Figures reproduced from reference [5]. 
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The EAL concept is used in a depth-profiling experiment, where the change in photon 
energy leads to a change in the kinetic energy of the ejected photoelectron, which variates 
the amount of signal coming from the surface as compared to the bulk. 


For aqueous samples, this universal curve is still poorly known. An important work [5] 
used the microjet technique from O1s photoemission intensity as a function of photon energy 
to determinate experimentally the curve for liquid water, figure 3.4b, which its minimum is 
estimated to be about 20 “A. This value is different to what was previously determined (5 - 
10ºA) [48, 49, 50]. 


Photo lonization Cross section 


The probability of a photoelectron being emitted from its electronic state is called 
photoionization cross section. The differential photoionization cross section (DPCS) 
provides the number of photoelectrons emitted in a given solid angle dO, in other words, this 
quantity gives information about how many photoelectrons are ejected in different directions 
concerning the incoming photon beam. DPCS depends on some factor such as the incident 
photon energy and polarization as well as which orbital in the target atom or molecule 
the ionization took place. In general, it is valid the proposition that for a photon with a 
wavelength greater than target atomic dimension the dipole approximation holds [51]. For 
an incident photon beam in the z direction with linear polarization vector in the x direction 
the differential cross section can be written as 


do o Ototal 


dO” 4r 


[1 + 8Po(cos6)] 
(3.2.2) 


Where 0,.,. is the total cross section and B is the angular distribution asymmetry 


total 
parameter [52] that depends on the photon energy and the symmetry of the ionized orbital, P, 
is the second-order Legendre polynomial. If ionization takes place in several core-orbitals or 
any atomic or molecular orbital depending on the direction we choose to collect the outgoing 
photoelectrons the relative photoelectrons intensity will change for different collection 
direction. There is, however, a special angle where the relative intensity of photoelectrons 
coming from different orbitals are the same as the relative intensity of photoelectrons if the 
each photoelectron intensity is integrated in all directions and then both relatives intensities 
are compered. This angle 0, in the dipole plane, is called the “magic angle” and is given by 
54.7 [53, 54]. This can be rationalized if it is realized that P,[cos(54.77] = O 


ELEMENT SENSITIVITY 


As previously stated, the photon energy determines which orbital is excited. X-ray 
photons excite core orbitals, and ultraviolet photons excite valence orbitals. Core orbitals, 
in first approximation, have atomic orbital feature and they are used for chemical analyses 
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of samples to identifying elements, due to the unique binding energy of core-electrons [55]. 
Valence orbitals are delocalized and participated in the chemical bonds; they have the 
molecular orbital feature and low binding energy. The water valence binding energy is used 
in this work as a reference to calibrate allthe CO, spectra and to monitor the stability and the 
conditions of the experiment, for this a photoemission spectra of liquid water is measured 
at the beginning of the experiment, and every time there is a sample exchange. Section 4.4 
presents more information about the energy calibration procedure. The binding energy of 
1b, liquid water is 11.16 eV [56], see figure 3.5. The peaks of liquid and gas phase present 
an energy chemical shift, although they represent the same energy level of the molecule, 
this shiftis a consequence of the different environment that the water molecules are in the 
liquid phase and the gas phase. The next example will present a complete explanation of 
the different chemical shift observed in a spectrum. 


For instance figure 3.6 below shows the spectrum of C1s ethanol 35%M. In this 
spectrum, there are two groups of lines, one for vapor structures and the other for liquid 
structures. The vapor structures have higher binding energies than liquid structures. This 
energy shift is the result of tree main contribution: electron polarization, surface dipoles and 
the water network of the H-bonding that generates changes in the orbital [57]. The liquid 
lines are broader than the vapor lines which is also an effect caused by the neighborhood 
molecules that generates different perturbations in each molecule that is being analyzed 
in the liquid sample. Thus the measured spectrum is an average of energy levels with 
small perturbations and therefore broader. The molecules in the gas phase are too far away 
from each other, due to this the electron binding energy does not suffer from neighborhood 
disturbances, and as a consequence, the spectrum line is narrower. From this figure, we 
can differentiate the carbon atoms that compose the ethanol molecule according to their 
binding energy. The hydroxyl group reduces the local electron density, resulting in a less 
effective shield on the carbon nucleus, so the core-electrons have a higher binding energy 
methyl C1s peak. 
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Figure 3.5: Photoemission spectra of liquid water valence measured at 360 eV photon energy 
that shows valence water orbitals and its vibrations. This spectrum is used to calibrate the 
spectra used in this work. The 1b, of liquid has a lower binding energy than gas phase. 
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Figure 3.6: A Cis spectrum of 35%M of ethanol in water. Figure presents signal from the 
liquid and vapour structures. More information in text. 


THE PHOTOEMISSION LINE PROFILE 
According to the uncertainty principle, there is a relation between energy E and time 
which is given by: 


SE(eV ). (fs) > 1 (3.4.1) 


Since there is no time operator in quantum mechanics, time, in this equation, 
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actually means lifetime T or, in another word, the time a system of interest takes to undergo 
substantial change. In the present case, this system is an atom or molecule undergoing 
photoionization. Equation 3.4.1 has a significant consequence for the present spectroscopic 
study. It settles a limit for the "intrinsic widths” of the core or valence spectral features. 
According to the uncertainty principle [58] peaks in the photoelectric spectrum cannot be 
infinitely narrow. In reality, the width of a peak, dE, which is a measure of how well defined 
the transition is determined regarding energy E - is ultimately defined by the lifetime of the 
states involved in the transition. For example, for an ionized core state, the lifetime of the 
ground state convoluted with that of the final ionic state will set a limit on the intrinsic core 
line width. Usually, the component coming from the ground state may be disregarded and 
the lifetime of the core ionic state will define the intrinsic width. It can be demonstrated that 
the dispersion in energy of the signal strength of a valence or core peak has a Lorentzian 
profile. In this case the only spectral broadening is caused by the uncertainty principle. In 
this thesis, the Lorentzian profile full width at half maximum (FWHM) of the Cis peak, is 
typically 0.10 eV [59]. FWHM is a measure of 5E(eV ) in the uncertainty principle formula 
3.4.1. According to that formula, the lifetime of this core-hole is about ten f s. 


So far it has been talked about the natural line width, the broadening resulted from 
the radiation process. Other broadening factors, present in an XPS peak, usually dominate 
its profile. Several factors can generate a Gaussian contribution to the profile of the line. 
Molecules, in a liquid, present a Maxwellian distribution of velocities. Those different speeds 
are given rise to the Doppler shift that will be broad the line [60]. This effect is temperature 
dependent. Another contribution to the Gaussian format is due to instrumental effects, mainly 
due to the monochromator slit that determines the photon energy bandwidth and the pass 
energy of the electron spectrometer [61]. The resolution of the experimental apparatus must 
be calibrated so that the width is given mainly by the width of what is being investigated. The 
significant contribution to the broadening comes from the chemical shift [62], as explained 
in the previous section 3.3. 


When all those effects are presented in a photoelectron spectrum, the resulting 
profile is the convolution of the Gaussian and the Lorentzian contribution which gives a 
Voigt profile [60] 


V(a) = E L(y)G(x — y)dy (3.4.2) 


oo 


Land G are the Lorentzian and Gaussian functions. The FWHM of a Voigt is given by 
an association of Gaussian and Lorentzian width. The FWHM of figure 3.6 for CH, is 0.573 
mevV, as said before the Lorentzian FWHM is 0.1 meV and the Gaussian FWHM is 0.517 
mevV. Therefore the Gaussian profile makes the major contribution on the FWHM of a Voigt. 
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Pos-collision interaction (PCI) is another factor that can affect the shape of XPS 
peak [63]. This interaction happens when core level photoionization is followed by Auger 
decay. Usually, the photoelectron kinetic energy is smaller than that of the Auger electron 
which will then pass the photoelectron. Under this situation the electrostatic interaction 
between the nucleus and departed photoelectron will be higher, resulting in a delaying effect 
depending on the photoelectron-ion distance. This distance will vary since the Auger decay 
is a stochastic effect. This stochastic delaying effect generates an asymmetry in the XPS 
profile toward higher binding energies. 
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CHAPTER 4 - MICROJET TECHNIQUE AND NMR 


In this chapter, it will be presented the experimental setup and the microjet technique. 
The synchrotron radiation and relevant specifications concerning beamlines where the 
experiments took place are also presented. The liquid jet thermodynamics is discussed. Itis 
also described the fitting procedure employed to the XPS data. Finally, it is reviewed some 
concepts belonging to the nuclear magnetic resonance spectroscopy (NMR) which was 
used to obtain information about the bulk CO, species concentration. 


SPECTROSCOPY ON LIQUID SAMPLES 


Itwas discussed before that the electron inelastic mean free path (EIMFP) is minimal 
under ambient pressure. Thus electron spectroscopy requires high-vacuum conditions. 
When dealing with solid or low-pressure gas samples, the EIMFP is long enough for the 
electrons to reach the detector without considerable losses. Performing XPS on the surface 
of volatile liquids, however, pose a real challenge due to the small EIMFP. For example, water 
has a vapor pressure of 23.3 mbar [20] at room temperature. Under vacuum conditions, 
evaporation from the liquid surface will keep a dense cloud of vapor molecules above the 
liquid surface increasing changes a photoelectron would be scattered inelastically thus 
reducing the EIMFP substantially. 


In 1973, Hans and Kai Siegbahn performed the first liquid phase spectroscopic 
experiments on liquid formamide. They demonstrated the complete separation between 
the liquid and gaseous features [64]. The experimental setup consisted of a rotating metal 
disc partially flooded into a pool of liquid sample. In this way, the liquid surface was always 
renewed. For this configuration, the sample should have a low vapor pressure or salt should 
be added to lower the vapor pressure of more volatile samples. 


Fifteen years later, M. Faubel proposed a new approach which successfully allowed 
the study of more volatile liquids by using a liquid micro-jet [21]. With this technique, the 
liquid is pushed by a glass nozzle forming a cylindrical jet of some micrometers that is 
injected into the evacuated experimental chamber. The jet leaves the glass nozzle with 
a high speed, which ensures that the liquid surface is continually renewed. There is a 
manageable increase in the pressure of the experimental chamber due to the liquid that is 
continuously injected into the cnamber and freeze further downstream using liquid nitrogen 
traps. This technique allows measurements of both the liquid phase and the gas phase. 


The cylindrical microjet can be divided into three parts, the liquid bulk and surface 
are situated in the innermost part region |, see figure 4.1b, region Il contains the vapor and 
lllis characterized by the free molecular flow. In the vapor region the ejected electrons may 
further lose their kinetic energy by inelastic collisions; therefore it is essential to minimize 
this area. One way to achieve this is to introduce a skimmer into this region. The effect of a 
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skimmer is to change the geometry from cylindrical to spherical, so the pressure will drop 
from 1/r to 1/12, which ensure a differential pump between the experimental chamber and 
the analyzer. 


(b) 


Figure 4.1: (a) Liquid microjet injected in a vacuum chamber with 18 um of diameters and 80 
m/s. (b) Regions of a cylindrical jet. In lis shown the liquid and surface region, in Ilis shown 
the vapor region and Ill is shown the free molecular flow region. 


To get a sufficient intensity of photoelectrons, the pressure of the liquid P and the 
distance d between the liquid surface and the entrance of the electron analyzer has a 
crucial role. For a cylindrical geometry of microjet, the relation between those parameters 
is approximate by 


Pd = PR. (InR 


o je 


-InR,) (4.1.1) 


spec 


where P, is the vapor density, R,, is the liquid jet radius and R . is the distance from 


e 


the center of the liquid jet to the skimmer entrance [65]. 


PHOTON SOURCE: SYNCHROTRON RADIATION 


When a charged particle near the speed of light undergoes a centripetal acceleration, 
the emitted radiation is called synchrotron radiation. Many synchrotrons around the world 
provide this type of beam which allows scientists to perform various experiments. Among 
them, Laboratrio Nacional de Luz Sncrotron (LNLS) in Brazil and MAX-lab in Sweden were 
used as a photon source for the present experiments. Synchrotron radiation has interesting 
properties which make it a great tool for the study of liquid micro-jets. Among them, we may 
mention high brightness, wide emission spectrum (from infra-red (10! eV) to hard x-ray (105 
eV)), a high degree of polarization, high collimation and a high stability under high vacuum 
conditions. There are other forms to obtain x-rays, such as, x-rays tubes or gas discharge 
lamps, but the photon energy of those devices is very restricted. 


In a synchrotron, the electrons are kept in storage rings made of vacuum tubes. 
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Normally the ring has a shape of a polygon formed by straight sections with dipole magnets 
in the edges changing the trajectory of the electrons, producing synchrotron radiation at 
specific points of the ring. Those bending magnets produce a broad bremsstrahlung. Other 
devices used to produce synchrotron radiation are insertion devices. They are placed in 
straight sections of the ring. For this thesis, the most important devices are the undulators. 
They consist of a periodic structure of magnetic dipoles with alternated polarity, forcing the 
electrons to make an oscillatory movement when they are passing through. A change in 
the gap between the upper and lower magnetic structure will change the distribution of the 
photons produced. 


In the beamline, the synchrotron radiation is brought to users. After leaving the 
undulator, a focusing mirror collects the non-monochromatic divergent radiation towards 
a monochromator. In which it is possible to choose a photon energy range depending on 
the beamiline and the type of experiment. The resolution is determined by the exit slit of 
the monochromator and the monochromator grid. Between the monochromator and the 
experimental chamber, there is a refocusing mirror that produces a final photon beam 
collimation. 


Beamlines and experimental setup 


The photoelectron spectroscopy studies presented in this thesis were performed at 
LNLS in Campinas, Brazil, at the beamline U11 PGM (Plane Grating Monochromator) [66] 
and MAXlab in Lund, Sweden, on the beamline 1411 [67]. Both beamlines are permanent 
multi-purpose end station, where experiments with gas, liquid, and solid samples is 
performed. The photon energy of interest for spectroscopy studies in liquids are in the range 
of 30 - 2000 eVs. At U11 PGM, the range of photons produced by the beamiline is 100 - 1000 
eVs, while at 1411 the range is of 50 - 1000 eVs. 


On both beamilines, there is a hemispherical electron energy analyzer Scienta 
R4000, but the slit in the entrance of the analyzer at U11 PGM used during the experiment 
was of 200ym and at 1411 the slit was of 50ym, the pass energy in both experiments was 
200 ev. All the experiments at 1411 were performed with a detection angle of 54.7º, the 
so-called “magic angle” as shown in section 3.2.1, concerning the polarization plane, the 
photon direction and the propagation direction of the liquid jet were perpendicular. 


Figure 4.2.1 shows the setup of the experimental station for liquid jet experiments. 
In the experimental chamber the stainless steel rod, where the source of the micro-jet 
is located, allows precise positioning of the jet concerning the entrance of the analyzer. 
As shown in figure 4.2.1, the synchrotron radiation intercepts the laminar flow region of 
the liquid jet, and this intersection should be below the skimmer which diameter is 1 mm. 
A Nitrogen trap collects and freezes the jet further downstream; this helps to maintain a 
manageable low pressure in the experimental chamber. Photoelectrons passing through 
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the skimmer toward the spectrometer find an area where electrostatic lenses focus that 
accelerate / retard them to the electron analyzer entrance slit, due to the high vacuum in 
this region the energy loss by inelastic collisions can be considered negligible. The pressure 
inside the analyzer cnhamber is between 103 to 10º mbar. 


A High-Pressure Liquid Chromatography (HPLC) pump leads the solutions through 
PEEK (Poly-Ether Ether Ketone) to a stainless steel tube and to the glass nozzle which 
had an orifice of 20ym. The pressure in the pump is adjusted so that the jet formed has a 
constant flow speed given by 50 m/s. This way it is possible to obtain a liquid jet in a vacuum, 
for more viscous liquids the pump pressure must be higher to maintain the constant flow. 
To switch between samples during experiments, it is used a system of beakers connected 
to a switching cross. 


Five samples were prepared with 0.5M of Na,CO, in demineralized water (Millipore 
Direct-Q) at pH 7.0, 8.8, 9.8, 10.8 and 12.4. The pH of the samples was adjusted with HCl 
or NaOH and the samples were filtered to remove particles that could clog the glass nozzle. 
After and before the spectroscopy measurement the pH of the samples was verified. A 
small amount of NaCl, 50 mM, was added to all samples to avoid charging effects caused 
by some sources among them two are worth mention: the photoelectric effect itself and 
electrokinetic charging. This last effect is connected to the charge build up between the 
glass nozzle small channel walls and the sample. Charging, if not suppressed will cause 
broadening and shifts in the photoelectron spectra. 


electron energy analyser 


my W€— Skimmer 
nozzle 


synchrotron radiation 


turbomolecular pump 
3x 300 I/s 


main chamber 


turbomolecular pump 
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Figure 4.2: Experimental setup of photoelectron spectroscopy on liquid microjet. 
Energy analyzer 


After the electrons pass through the skimmer, they will reach the electrostatic lenses, 
which focus the electrons in the entrance slit of the hemispherical analyzer. The electrostatic 
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fields inside the analyzer select the electrons with a given kinetic energy that will reach the 
middle of the micro-channel plate detector, this energy is called the pass energy. Depending 
on the pass energy the electron-lenses may accelerate or retard the electrons, according to 
their initial energy. In average they will have the same energy as the assigned pass energy 
at the analyzer entrance slit. The spectrum is recorded by scanning the voltage applied to 
the electron lenses. 


The slit length w, the pass energy E, and the analyzer diameter 2R set the resolution 
as 


Ep 
AE xu- E (4.2.1) 


In practice the pass energy and the slit size are the parameters that can be changed, 
itis important to find a suitable combination of them to reach the best possible resolution. 
There are two multi-channel plates (MCPs) inside the detector. Their function is to multiply 
the electrons by a factor of up to 10”. After multiplication, these electrons hit a phosphor 
screen. The impact results in bright spots that are recorded by a charge-coupled device 
(CCD camera) located outside the detection cnamber. The CCD camera gives information 
about kinetic energy and spatial/angular distributions. 


SURFACE EQUILIBRIUM 


Several experiments were performed to elucidate the surface thermodynamic 
equilibrium of a micro-jet at the measuring position [62, 68]. Water has a relatively high 
vapor pressure and, in a vacuum, the evaporation process is intensified, therefore at the 
measurement position, the liquid surface must be well behaved. 


detector junction lens 


boxes 


(a) (b) 


Figure 4.3: (a) Electron analyzer. The electrostatic lenses and the detector are depicted [6]. 
(b) Details of the hemispherical analyzer. w is the entrance slit, R, is the radius of the internal 
sphere, R, is the radius of the external sphere and R is the average path of the electrons. 
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Another import point is the diffusion time of species within the microjet. In principle, 
one may question if the time to reach a steady state is longer than the time the jet takes 
to reach the ionization region. The issue about the evaporation process is answered by 
theoretical assumptions. Molecules may condense or evaporate from a liquid surface at a 
rate per unit area j given by: 


. 1 SkT me Pgas 
7= (Vgas) Ngas = 4 ( ) ( e) (4.3.1) 


Megas 


where Mada? is the average speed of gas particles and Nas is the number of gas 


s' 


particles, on the right side equation the m,., and p,., are the mass and the partial pressure 


as 
of the gas respectively. K is the Boltzmann si and T is the temperature [62]. 

The liquid water monolayer completely evaporates after approximately 100 ns at 8 
mbar and 0ºC. The liquid molecules collide with each other on a time scale of 1ns, these 
collisions lead to thermodynamic equilibrium within 10ns [62]. The liquid-liquid collisional 
relaxation is at least 4 orders of magnitude faster than the collision events of the gas phase. 
A jet with a velocity of 100m/s reaches the electron beam after 10y's. So at the moment of 
measurement the liquid surface of the microjet is already in thermodynamic equilibrium [69]. 


Experiments with two different setups and different pressures but with the same 
sample were performed [48]. No difference was observed between the two experiments, 
therefore the low pressure surrounding the micro-jet must be regarded as not being a 


relevant factor. 


The diffusion kinetics of a solution determines the formation time of an adsorbate 
monolayer. Aqueous solutions at the concentration limit were investigated and it was not 
possible to observe variations in the surface concentration as a function of the distance 
between the nozzle exit and the observation region [62]. 


CURVE FITTING 


The data treatment and including the fitting procedure are important tools for the 
correct interpretation of a spectrum. For such activities, a knowledge of physics behind 
the process is necessary to obtain meaningful information about the spectrum under 
consideration. In this thesis, the SPANCF (Spectrum Analysis Curve Fitting) package 
was used to perform all the data analysis [70]. The package is a plugin to IGOR Pro by 
Wavematrics [71]. 


Before starting the fitting itself, it is necessary to normalize the spectrum. In this 
way, intensities between different runs is compared. Normalization is made concerning the 


photon beam intensity, which it decreases over time and number of sweeps. The ring current 


Chapter 4 - Microjet technique and NMR ES 


is proportional to the photon beam intensity, so sometimes it is sufficient to use the ring 


current information recorded during the experiment. The normalized intensity 1, is given by 


Tx = 1000 - 
N ae (4.4.1) 


where the n is the number of sweeps, c, is the ring current and the multiplication 
factor 1000 is used to have high enough values, to avoid problems with truncations. 


The energy axis needs to be calibrated before starting the fitting. Thus small energy 
shifts may be present due to variations in the experimental conditions. The calibration 
is typically performed using a known peak, such as the 1b, transition in the liquid water 
whose binding energy is 11.16 eV [56]. For the calibration procedure and also intensity 
normalization a water valence spectrum is recorded before and after a measurement. By 
doing so any small change in energy as well as intensity is monitored and corrected for. 


For the fitting itself it is necessary to define initially the parameters that characterize 
the spectra. The first step to fitis to choose the line-shape, the chosen one for all the peaks 
of this thesis is the Voigt, which is a convolution between a Lorentzian and a Gaussian line 
shape. There is no PCI asymmetry in the system, so this parameter is fixed to be zero to all 
specitra. The SPANCF package allows to set parameters during the fit or to let then variate 
free. Itis also possible to do links between parameters from different spectra. Those links 
are essential to extract relevant data from the fitting. There are three peaks in each spectrum 
of this thesis, which represents carbon dioxide (CO), bicarbonate (HCO, and carbonate 
(CO?). The energy from CO, peak was fixed to be the same as the pH 7.0 spectrum for all 
spectra. This peak was chosen because the CO, only appears on this spectrum and even 
thought it is not significantly present in the other samples it is important to maintain this 
peak for the fitting of the other spectra. The energy of HCO," and CO ?-was linked in a way 
that the energy difference between both peaks is the same as that of pH 9.8 for all spectra. 
This constancy in energy is plausible because it is only the concentration of the species 
which varies from sample to sample, so the binding energy should be the same. This link is 
relevant because some spectra presents only one of the species, the signal from the other 
one is below the background and therefore fixing the binding energy difference allows that 
the peak of lower intensity is properly adjusted. The 9.8 pH spectra was chosen because 
the intensity of both species is meaningful.The intensity of all peaks is free to variate. The 
Lorentzian FWHM (full width at half maximum) is fixed to 100 meV [59], and the Gaussian 
FWHM of each species was free but made to be the same for all spectra. 


To compare intensities between different samples the experimental setup has to be 
stable during the measurements. In sections 3.2 and 3.2.1 it was said that the intensity of 
the photoelectron signal is proportional to the number and the cross-section of the species 
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investigated in a specific volume. From the ratio of the intensities from different samples it 
is possible to obtain information about abundance and propensity of species to the surface. 


SURFACE ENRICHMENT MODEL 


In section 2.3 it was discussed the Hendersen-Hasselbalch equation that describes 
the behavior of acid-base reactions in bulk. The behavior of those reactions on the surface 
must be similar to that described by equation 2.3.9, however, this model needs to be 
improved to describe correctly the species variation on the surface. Due to the presence 
of the surface, some molecules are more or less prone to stay at the surface, which may 
increase or decrease the concentration of those species differently ftom what is found in 
bulk. This surface enrichment factor g needs to be added to the Hendersen-Hasselbalch 
equation as proposed in [72]. 


For this, the photoelectrons intensities for both species is rewritten as 


[HCO; = x.[HCO3] (4.5.1) 


Ihcoz 


[cos] = 9,02-*-[005] (4.5.2) 


So the intensities are given by the enrichment factor g, for each species, by an 
experimental factor x that depends on the experimental setup and by the concentration of 
each specie. This experimental factor may be ignored if the intensities are taken at the same 
experimental setup. The acid fraction in the surface is writing as 


JHCO; HCO] 1 


JHcos [COS] 4 Icoz-[C05] 1+ [coz] (4.5.3) 


gHCO, ] 


where the ratio between the enrichment factor g is defined as g = 9ucoz 1902 . Using 
3 
equation 2.3.6 in the equation above 


1 


Rs = 102H-pKa-log(g) +1 (4.5.4) 


This expression is similar to 2.3.9, but the g factor shifts the curve depending on the 
concentration of each species in the surface. For instance, the curve is shifted to higher pH 
if the acid is more enriched; otherwise, the curve is shifted to lower pH if the base is more 
enriched in the surface. An apparent pK, can be defined from equation 4.5.4 as 


PK a = PK, + log(g) (4.5.5) 


a,app 
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The molecules studied in this work are charged and the surface behavior as an 
interface that separates two dielectric mediums. The punctual charges in one of the 
mediums generate an image charge of the contrary signal on the other medium. In this 
way, the interaction between the charge and your image charge difficult the creation of a 
monolayer of ions on the surface. Therefore it must be added a term lin the proposed model 
that corrects the fit when the solution is saturated with a certain species. Therefore equation 
4.5.6 presents the acid fraction as a function of solution pH for the surface. 


4 
= 107H-pKa-log(g) +1 (4.5.6) 


Rs 


From the fitting performed using this equation, it can be done a comparison with the 
results obtained by the Hendersen-Hasselbalch and observe the similarities and differences 
between both. 


LANGMUIR MODEL 


If the studied species have a high propensity to populate the surface as compared 
to the bulk, its corresponding photoelectron signal, having a kinetic energy amounting to 
about 100 eVs, will have a higher intensity than that from species depleted from the surface. 
In this section, we will explain the formulation used in this thesis to estimated the surface 


concentration. 


The Langmuir model describes systems performing ideal chemical adsorption, and 
it is derivative from basics models containing kinetics, thermodynamics, and statistical 
mechanics. The model was established to describe the adsorption of gas molecules in solid 
surfaces. In the case of liquids, the molecules are considered to be adsorbed to the surface 
from the bulk when it has a high propensity to the stay on the surface. There are three 
assumptions to apply this model 


*. The adsorption occurs in specific sites, and they are all identical; 
*. The adsorption energy is independent of how many sites are occupied; 


*- Qnly one adsorbed molecule occupy each site, and once all sites are occupied 
the adsorption ends, so we have the formation of a monolayer. 
There are just a few systems that have this kind of behavior, but despite all those 
limitations the Langmuir model is a perfect first approximation and is widely used to fit 
experimental data from XPS in liquids [73]. 


The total signal is composed by a surface part |; and a bulk contribution |,. Itturns out 
to be more convenient to use the bulk molar fraction x, instead of |, from the experimental 
point of view. In this way the Langmuir isotherm is 
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Loma “Xb 


Is = : (Saga) 
Where the | is the maximum photoelectron signal expected from the surface 


s,max 


molecules when a monolayer is formed and AG,,, is the Gibbs free energy. These two 
parameters are determinate with the fitting procedure. The temperature T during the 
experiment was 283.15 K and the universal gas constant used is R = 8.31446 J/(K mol). 


Usually the Langmuir isotherm is built in a way that the signal intensity variates as 
a function of bulk molar fraction at fix pH. In the case of the reaction studied in this thesis, 
changing the pH of the solution alters the chemical equilibrium. In this way the Langmuir 
isotherm is approached in another way, using pH change to generate the variation in the 
bulk and thus to observe the signal intensity of photoelectrons. 


A full description is presented in section 5.4, where the data are also presented. 
This new approach, proposed in this thesis, could be applied for others molecules that 
deprotonate at different levels of acidity. 


NUCLEAR MAGNETIC RESONANCE 


Nuclear Magnetic Resonance (NMR) is a well-known technique so that this section 
will be a short review of the principals and important physical aspects of it. NMR allows 
obtaining information about the concentration of different species with varying degree of 
protonation in a solutions bulk which, in our case, is very useful to interpret the results 
obtained by XPS. NMR involves the transition between states characterized by the different 
orientation of nuclear spins subjected to magnetic fields. The transitions are induced by 
electromagnetic radiation. 

The nucleus has intrinsic angular momentum | and a magnetic dipole associated q. 
The element of interest for this thesis is carbon that has a null nuclear spin, but naturally, 
there is 1.109% occurrence of carbon 13 that has spin 1/2, these were the carbons used to 


perform the experiment. 


1 1 y 
= (; + 1) Sha 8=(1+1) 18 (471) 
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a = 50 «B=-—28 (4.7.2) 


The magnetic moment in terms of lis given by 


q , 
H= Nom! = gnBnl = nt (4.7.3) 


gy is the nuclear g factor, B, is the nuclear magneton, m, is the nuclear mass and y 
is the magnetogyric ratio. 


In the presence of a magnetic field, the nucleus is either aligned or anti-aligned. The 
difference in energy between an aligned nucleus with the magnetic field in the z direction 
and one against it is of AE =hyB.. If this nucleus is irradiated with electromagnetic radiation 
at the same frequency of AE, that is AE = hv, the radiation will cause transitions between the 
states aligned and anti-aligned. The frequency associated with transition of the aligned state 
to the anti-aligned is directly proportional to the magnetic field, given by 


v= toe (4.7.4) 

The experiment can be performed by fixing the magnetic field and varying the radiation 
frequency. A signal will be detected when an absorption occurs. Another possibility is to 
fix the radiation frequency and scan the magnetic field. Both methods are equivalent and 


present the same results. 


The sample is placed between poles of a strong electromagnet, as in figure 4.4. The 
magnetic field is then varied by changing the current around the electromagnet. The sample 
is irradiated by an electromagnetic wave for which its frequency is scanned by 4.7.4 the 
absorbed radiation by the sample is detected and recorded. 


Depending upon the chemical environment, neighbor molecules may affect the 
nuclear energy which generates a chemical shit. This chemical shift has the same meaning 
as the one explained in section 3.3. In this case, the electrons of the molecule, in the 
presence of a magnetic field, generate an additional magnetic field that shields the nucleus. 
This shield depends on the electronic and/or chemical environment. Therefore depending 
on the environment the electromagnetic frequency in which the spin transition occurs will 
be slightly different 


27v 
Y(l-o) 


where B, is the magnetic field that the nucleus effectively exposed to and o is the 


Bo = (4.7.5) 


shield constant. This is the way that the chemical shift is presented in NMR [74]. 
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Figure 4.4: Schematic diagram of a magnetic resonance spectrometer. 


This technique of continous-wave (CW) is not efficient as the Fourier analysis 
technique in which all the nuclei in a molecule are excited at the same time with more 
than one frequency, when the molecules relax to equilibrium state the frequencies emitted 
are analyzed. The samples are exposed to very short and strong pulse of a given carrier 
frequency that contains a range of frequencies along one axis all the nuclei of interesting 
are excited. The magnetization vector in this case is the sum of all nuclei excited. The signal 
is an overlapping of all the nuclei relaxing, they are collected and followed by a Fourier 
transform which converts the complex time domain signal emitted into the frequency domain 
spectrum. From this method it is possible to obtain a full spectrum faster than CW. As a 
consequence of the relaxation mechanism, that is first order process, the signal emitted 
decays exponentially, which is known as free induction decay (FID). For samples that 
presents low abundance isotopes as 'ºC the FID's signal is weak. In this case the FID's 
signal may be stored until it is enough to provide a spectrum. 

The NMR experiment was performed at the Brazilian Biosciences National Laboratory 
LNBio in Campinas, Brazil [75]. A small part of the samples was transported from the XPS 


experiment in Sweden to Brazil so both experiments are performed in the same samples. 


Curve fitting 


Previous work demonstrated that it is possible to determine the pH and solute 
concentrations from aqueous solutions of CO, by "ºC NMR [76]. From the NMR spectra, 
it is possible to obtain the ratio between concentrations of HCO, and CO,? species. The 
chemical shift depends on the relative concentration of those two species. At low pH, the 
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solution has 100% of protonated species (HCO), so the chemical shift is Vucoz . Athigh pH 


the solution has 100% of deprotonated species (CO,?), so the chemical shiftis V.o2- At the 


cos” 


pH where both species are present in solution, the chemical shift is given by v 


V=Vico;” Xnco; + VYooz- Xcog- (4.7.6) 


where Xsco; ANd Xcog- are the molar fraction of the two species at the same pH. 
The molar fraction of a specific pH of the deprotonated specie is determined by using the 
chemical shift as 


chemical shift at low pH — observed chemical shift 


Xooj” = chemical shift at low pH — chemical shift at high pH (tao 


Then using that Xuco; + Xcog- = 1itis possible to evaluate the molar fractions of the 
protonated specie. After the determination of both molar fraction the pK, is obtained using 
a version of 2.3.6 


X - 
DK, = pH +log HCO; 


Xcoz- 


(4.7.8) 
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CHAPTER 5 - RESULTS AND DISCUSSION 


This work presents a study of the molecular interface of aqueous solutions relevant 
to the environment. More specifically, the concentration of two CO, protonated species at 
the molecular surface dissolved in liquid water is investigated. 


In chapter 2 it was pointed out that the distribution of solutes in bulk is homogeneous, 
on the surface however this distribution may differ substantially depending on some factors 
which it is expected to elucidate some of them for the present solution. For example, 
biomolecules pursuing a more extended carbon backbone tend to concentrate more on the 
surface [77]. Solute molecules residing in the aqueous molecular surface may generate, 
besides others changes, alterations in the hydrogen bonding network, which is an active 
subject of study in the community [78]. 


Among various characteristics of x-ray spectroscopy, those that are the most 
important to the present study are the element specificity, the surface sensitivity when the 
photoelectrons energy is correctly chosen and the chemical state sensitivity, as explained in 
chapter 3. It will be presented in this chapter the results from XPS experiments performed to 
comprehend the distribution of protonated molecules in the aqueous interface. Two sets of 
XPS measurements were performed using a microjet, however, in one of the setup defective 
power supplies prevented us from obtaining the desired high-quality spectra needed to 
perform an unambiguous study. Thus, the data related to this experiment will not be 
presented in this thesis to avoid confusion. The XPS measurements were complemented by 
data from another technique namely nuclear magnetic resonance spectroscopy (NMR). In 
it, the same samples were investigated. Since this technique is only bulk sensitive and also 
provide a way to probe the distribution of protonated species in that region of the solution 
it was used to eliminate some spurious effects to explain our observations. The results of 
both experimental techniques are compared as a function of pH, and the observed changes 
are discussed. 


Solutions that contain more than one solute are widespread in nature. For example, 
the ocean water contains ions such as sodium (Nam), chlorine (CI), magnesium (Mg?), 
calcium (Ca?), potassium (K*) and sulphate (SO%) in the form of dissolved salts. Apart from 
these, there are other elements like carbon (C), nitrogen (N ), phosphor (P ), carbonate 
(CO,2) and oxygen (O,) which are essential to the growth of marine living beings. One must 
bear in mind however that the experimental investigations of these solutions with more 
than one solute is quite challenging. One reason for this is that is quite hard to distinguish 
between interactions among solute themselves. In this study three solutes will be dealt with: 
bicarbonate HCO,”, carbonate COZ- and the carbonic acid H,CO, that in its majority will be 
present in the gaseous form of CO,. The species studied are charged, so it is expected that 
they tend to reside in bulk of the solution [40]. 
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THE SURFACE 


As we pointed out in the introduction the CO, gas present in the atmosphere may 
absorb at the surface of oceans and aqueous aerosols in the atmosphere. Once inside the 
liquid, it reacts with water molecule forming carbonic acid, and depending on the solutions 
acidity, this acid may be deprotonated forming bicarbonate or carbonate with pK, 3.6 and 
10.3, respectively, both are liquid species. The fact that both species are liquid usually is 
pointed out as one of the main reasons why CO, can be dissolved in aqueous solutions 
to a relatively high concentration. Depending on the amount of CO,/H,CO, dissolved in 
a solution, the acidity will change and, as a consequence, the distribution of H,CO, 
deprotonated species in the solution will also change. For instance, at pH higher than 10.4 
the majority of the molecules in the solution are doubly charged. A previous XPS study of 
aqueous amine-based solutions and the process of CO, capture using a liquid micro-jet 
showed that the compound formed by the CO,, water, and the amine tends to go to the bulk 
[79]. It does remain however to be elucidated how the deprotonated species discussed here 
behave in pure water. 


The hydrogen bonding network is altered by the presence of a solute in the solution 
[78], this fact may lead to a decrease in the surface tension compared to pure water. This 
phenomenon explains in part the propensity of a compound. Compounds reside closer to 
the surface alter the solute distribution in the interface. The vacuum interface can be seen 
as a hydrophobic surface highlighting hydrophobic interactions, once hydrogen bonding 
cannot be formed in this region. The molecule orientation at the surface will not be treated 
in this thesis, but it is also an important factor, once different orientations result in specific 
adsorption behavior to minimize the interactions between water and hydrophobic molecules. 
Some others interactions also increase the propensity to the surface of some molecules 
such as van der Waals interactions. 


The surface propensity 


As discussed previously, two sets of XPS data were acquired. However only one 
turned out to be usable for reasons we already discussed. The measurement performed at 
the Swedish synchrotron facility, MAX-Il in Lund at beamline 1411 [67] turned out to have 
excellent stability regarding intensity kinetic energy of the ejected photoelectrons. Both of 
these stated properties were checked by running water valance spectra in between each 
new sample of interested. In these valence spectra, the 1b, peak belonging to the liquid 
water was normalized, and the area and position of the peak were compared through the 
series of measurements. The energy stability of the valence 1b, peak was beiter than 50 
meVs, and total intensity showed variations smaller than about 5 percent. A few sources of 
stability could be listed for this beamline: the photon beam size is about 1 mm in the vertical 
direction at the focusing point while the microjet height was smaller than 25 micrometers, 
alignment of the microjet to the photon beam and the small analyzer nozzle was made very 
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stable by robustness of the alignment system of the analyzer cnamber. 


The sample preparation was explained in section 4.2.1. The photon energy was 
360 eVs which provided Cis photoelectrons with about 70 eVs, so the photoelectrons 
are emitted mainly from the first layer. The curve fitting for the data acquired with XPS is 
presented in section 4.4, here it will be presented the fitted spectra. Figure 5.1 shows all the 
spectra obtained for each of the five samples prepared. 
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Figure 5.1: C1s PE spectrum of 0.5 aqueous Na,CO, of five samples at different pHs. 


There are up to three Cis peaks in each spectrum: CO, gaseous (not present in 
all spectra), bicarbonate HCO,” and carbonate CO, each peak was fitted using Voight 
profiles colored as red, blue and green respectively. For all the five spectra, the binding 
energy belonging to each species tuned out to be the same as expected. Obviously, 
the concentration varied according to the sample pH. The binding energy of the K shell 
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of an isolated carbon atom is 284.2 eV [55], so the binding energies of the carbons in 
the molecules are shifted to higher binding energies, because when carbon performs a 
chemical bond to an oxygen atom. The oxygen electronegativity compared to the carbons 
attracts the carbon valence electrons. Therefore the carbon atoms become more positively 
charged thus attracting to a larger degree the outgoing photoelectron, therefore, their kinetic 
energy will be reduced compared to an isolated atom. For a fixed photon beam energy this 
means the binding energy of the C1s core lines will be increased. The binding energies 
found as a result of the fitting are given in table 5.1 with their corresponding FWHM. With 
the deprotonation, the binding energy of the carbon K shell is lower due to a build-up of the 
negative charge in the carbon atoms. This chemical shift is enough to make the identification 
of each peak in the spectrum. Therefore CO,> is the most negatively charged species and 
correspondingly posses the lowest binding energy. Next towards higher binding energy the 
HCO, The chemical shift between the two species is as much as about 1 eV, a typical value 
found for this type of change in the protonation state. 


Binding energy (eV) FWHM (meV) 
Co, 296.41 213.3 
HCO,- 292.85 1047.1 
co 291.80 1170.9 


Table 5.1: The binding energy and the FWHM of carbon dioxide, carbonate and bicarbonate 
get after the fitting procedure. These values correspond to the figure 5.1. 


The pH 7.0 spectrum is the only one with CO,. The red peak is identified as carbon 
dioxide and not as carbonic acid H,CO, because given the low pressure and temperature 
the majority of CO, is not converted into carbonic acid. Thus the signal contribution is mostly 
given by the gas carbon dioxide. Another evidence that this peak belongs to CO, gas is that 
itis substantially narrower than the others according to table 5.1. 


The peak areas are proportional to the concentration of the respective molecules 
at the surface given that the photoelectrons with about 70 eVs are representative of the 
surface region. Figure 5.2 shows the calculated acid fraction of each spectrum plotted as 
a function of pH. In this figure the experimental data is shown with a fitting curve given by 
the equation 4.5.4 based on the enrichment model. The black curve shows the acid fraction 
expected for the bulk using the Hendersen-Hasselbalch equation 2.3.9. The curves are 
shifted from each other as it is expected. 


In section 4.5 it was explained the reason for the | parameter in the curve fitting 
procedure. Figure 5.2 shows that there is no formation of a bicarbonate monolayer on the 
surface, even for low pH, when the concentration is saturated. In a similar fashion at higher 
pH, the curve does not tend to zero, which indicates that there is some remaining amount 


Chapter 5 - Results and discussion 


of bicarbonate, differently from what is predicted for bulk. Another information obtained from 
the data analysis is the surface enrichment factor. From the value of the enrichment factor, 
one can conclude that the doubly deprotonated species is less prone to stay closer to the 
surface than the singly deprotonated. This changes between doubly and singly deprotonated 
species can be rationalized by taking in to account that the charged species close to the 
surface form a charge image which brings as a consequence a net force pushing the charged 
particle towards the bulk [40]. The higher the charge the bigger the force. This general and 
simplified theory must be tested experimentally since it was found that lager and polarizable 
inorganic anions as Br and | enhanced at the surface [80]. 


The pk, value obtained by the fitting procedure is different from the one given in the 
literature since this value depends on several factors. In Section 2.4.1 it was explored the 


e C1s0.5MXPS HCO; 
—— pKa bulk : 10.3 
— Fit 
| = 0.92 + 0.04 
g = 1.0684 + 0.0001 
pKa = 9.8599 + 0.0004 


Acid fraction 


9.0 9.5 10.0 10.5 11.0 11.5 12.0 
pH 


Figure 5.2: Acid fraction as a function of pH. In red are the experimental data and the fit for 
the surface enrichment and in black it is shown the acid fraction expected for the bulk. 


pK, dependency with temperature, which is given by the van't Hoff equation. This 
equation presents a relation between the constant dissociation and the enthalpy of a 
reaction. The signal of the enthalpy gives the amount of heat absorbed or emitted by the 
transformation of reactants into products at a given temperature and pressure. So, for a 
positive enthalpy, the reaction absorbs heat, which means it is an endothermic reaction, for 
a negative enthalpy the reaction releases heat, which means it is an exothermic reaction. 
The deprotonation reaction studied in this work is an endothermic reaction and therefore for 
this type of reaction the value of the dissociation constant increases when the temperature 
decreases. The pK, value given by the literature usually refers to room temperature. In this 
case, the microjet is at a much lower temperature, about 10ºC. Taking into account only the 
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jet temperature the pK, value should be higher than the value given by the literature for bulk, 
but this is not what was observed with the help of the fitting procedure. As a conclusion, the 
temperature is not the determining factor for the observed pK, shift. 


The carbonate fit is performed with the basic fraction defined as 


R,=1-R, (5.1.1) 

Using this simple formula another fitting curve was obtained for the basic fraction as 
a function of the pH. The final result is presented in figure 5.3. From this figure, it is possible 
to observe at first sight that the species present a behavior similar to what is expected for 
the bulk, figure 2.5. The discrepancies will be more evident when comparing these data with 
those found for the bulk. 


The data referring to pH 7.0 was not used in the fitting because in that pH there is 
the presence of a third specie, gaseous CO,. The used model for the fitting only takes into 
account two species, therefore the presence of this third specie is not supported, and the 
curve obtained will be incorrect. Furthermore, the amount of Co, presented in the solution, 
is not correctly evaluated by the area of the corresponding peak at pH 7.0 since this peak 
comes possibly in totality from the CO, evaporated from the solution and overlapping the 
photon beam. Although this gas phase CO, possibly bears a proportionality ratio to the 
dissolved CO, several, difficult to control, factors influence this proportionality such as the 
way the sample is prepared, the actual temperature of the micro-jet, the time between the 
preparation of the sample and the measurement. All this points out to the unreliability of 
taking the CO, as the amount of dissolved CO,. The other two species being liquid are 
immune to all these uncertainties. 
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Figure 5.3: Bicarbonate (HCO;) and carbonate COZ concentration variation on the liquid 
surface as a function of pH. 


Chapter 5 - Results and discussion 


THE BULK 


Photoelectron spectroscopy is a powerful tool to study the species distribution on the 
surface, but this information must be compared to the distribution of those same species 
in bulk of the solution. For this purpose, we used other powerful technique, the nuclear 
magnetic resonance, on the same five samples of the previous experiment with 5% of 


deuterium dissolved. 


Figure 5.6 shows the NMR spectrum on the same five samples of the XPS and 
figure 5.5 shows the obtained chemical shift for carbon 13 using NMR for CO,, HCO, and 
CO, as a function of pH. The deprotonation reaction equilibrium between carbonate and 
bicarbonate is fast, so the carbon 13 spectra have only one singlet with an average chemical 
shift belonging to both species, this was not the case for the XPS measurement because 
those measurements represent an instantaneous picture of the situation at the aqueous 
solution. The photoelectron leaves the sample in a time smaller than femtoseconds. The 
NMR chemical shift depends on the concentration ratio of the two species. At lower (higher) 
pH, the signal is dominated by bicarbonate (carbonate). As the pH increased the chemical 
shift also increased from 162.3 ppm ("pure” HCO, to 168.1 ppm ("pure” CO,?). The fitting 
process shown in figure 5.5 was done according to what was described in section 4.7.1. The 
pK, value obtained using this process was 10.13, that is, the pH value where the derivative 
changes signal. 


One can readily notice that the pK, value found using the NMR technique is higher 
than the one obtained using the XPS data. 


It is also described in section 4.7.1 the process to extract the molar fraction of 
carbonate and bicarbonate from the chemical shift as figure 5.6 shows. This figure presents 
the variation in concentration as a function of pH of both species. The same behavior 
obtained here was found in those reported in the literature, figure 2.5. The increase of 
the pH leads to a decrease in bicarbonate concentration and an increase in carbonate 
concentration. The 10.13 pK, value is also exiracted from this fitting by the Hendersson- 
Hasselbalch equation and corresponds to the pH where the intersection between the two 
fitted lines, in other words, the pH where the concentration of both species is the same. 
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Figure 5.4: 13C NMR spectrum of 0.5 aqueous Na,CO, of five samples ar different pHs. 
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Figure 5.6: Concentration of bicarbonate and carbonate as a function of the pH in bulk. The 
triangles are the experimental data, and the continuous line is the fitting obtained from the 
Hendersson-Hasselbalch equation. 


COMPARISON 

In sections, 5.1 and 5.2 the concentration and variation of bicarbonate and carbonate 
species as a function of the solution acidity for surface and bulk regions were found. The last 
step left is to compare both and evaluate possible differences between them. 

Figure 5.7 presents data and lines belonging models discussed before and presented 
in figures and 5.6. Each color represents a molecule and an experiment type as it is explicit 
in the figure legend. 

The concentration of both species differs between bulk and surface regions, 
especially for higher pH. In this region there is a saturation, almost all the bicarbonate was 
transformed in carbonate; however, there is still some amount remaining of the product. 


Chapter 5 - Results and discussion DE 


The saturation is not complete because both species considered are charged so that they 
are repelled from the surface by the image charge created on the other side of the surface, 
forming a dipole. This image charge effect is more apparent at increased concentrations, 
and itis revealed in the value of the | parameter of the surface enricnment model. 


Another essential point to be noted is the pK, variation: both experiments were 
performed at different temperatures, NMR was performed at room temperature while the 
microjet temperature is about 283K. In this way, since the pK, parameter is temperature 
dependent, a change was expected. As discussed in section 2.4.1, the studied reaction 
is endothermic, and as the temperature decreases the pk, increases, this is not what is 
observed in figure 5.7. After this brief recap of what has been said in previous sections, it 
is important to mention that the adopted model takes into account only two species. For 
lower pH, a third species will affect the reaction. The presence of this third specie will cause 
a reduction of the HCO, species since this species is converted to the third species and 
forcing an apparent reduction of the observed pK,. A situation completely differently is found 
in the was not detected. Future molecular dynamics could help clarify this issue. 
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Figure 5.7: Comparison between the variation of the concentration of bicarbonate and 
carbonate as a function of pH. The balls are the XPS data, and the triangles are the XPS 
data. 


From the parameters obtained from the fitting procedure, it is observed that both 
species are prone to go to the bulk. The amount of CO, revealed by the XPS technique is 
not reliable because this signal comes from the gas distribution around the liquid micro- 
jet. Many parameters will affect amount evaporated CO, such as time elapsed between 
the preparation of the sample and the measurement, whether the sample was exposed to 
some degree of mechanical vibration and the XPS signal intensity variates strongly with 
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the alignment, but the other two species are liquids and free from these uncertainties. So it 
makes no sense to adjust the curve in the region for low pH. 


LANGMUIR ISOTHERM 


The Langmuir isotherm describes the adsorption of atoms or molecules into the 
surface. In the literature, this curve is constructed by taking the XPS signal as a function 
of solute bulk concentration. In other words, the experiment is performed by increasing the 
concentration of solute and monitoring its XPS signal (C1s peak area) as a function of its 
bulk Molar percentage. 


In this work it was observed that the change in the solutions pH shifts the reaction 
chemical equilibrium making the products to be favored in basic medium and the reagents 
to be favored in acid medium, thus varying the species concentration. Therefore, instead of 
making several samples with different concentrations and the same pH, it could make several 
samples with the same initial concentration but different pH that shifts the equilibrium. In 
this way the Langmuir curve is approached differently, using pH to vary the concentration 
and thus observe the intensity of the XPS signal. In this case, the initial concentration is the 
same for all samples. Thus it is difficult to determine the bulk concentration theoretically, but 
it can be experimentally evaluated using NMR since this technique is bulk sensitive. At the 
end of this procedure, the Langmuir isotherm will be fitted to the XPS signal as a function of 
bulk concentration as is usually done in literature; the novelty is the way the concentration 
changes. 


The results for carbonate and bicarbonate is shown in figures 5.8 and 5.9 
respectively. For carbonate equation, 4.6.1 fits very well, and it is possible to make a good 
work determining the Gibbs free energy to be - 3.7 + 0.9 kJ/mol, see figure 5.8. | max is the 
other parameter to be determinated by the fitting, which is the photoelectron signal observed 
when the surface is completed saturated. This Gibbs free energy may be compared with 
that obtained recently for Ethanol which was -10.8 kJ/mol, which forms a monolayer at 2 mol 
percentage [78]. Another interesting value to analyze is the obtained |. 1, = 1.9 which is ten 
times larger than the species concentration at pH 12.4 meaning that in reality a monolayer 
will never be formed. This last conclusion is in line with the fact that this species avoids the 


surface. 

Bicarbonate presents the reverse behavior of the previous case. However, equation 
4.6.1 does not fit the data very well, and the Gibbs energy is nearly zero, figure 5.9. This 
may be explained by the fact that bicarbonate participates in two chemical reactions, see 
figure 2.6, which is not taken into account by the model. 

This proposed way to apply the Langmuir model may be a useful tool to study others 
species undergoing deprotonation process such as amino acids, figure 3.2. 
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Figure 5.8: The points represent the photoelectron signal intensities |, of carbonate as a 
function of your bulk molar fraction x, - The variation of concentration is given by the 
different pH. 
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Figure 5.9: Langmuir isotherm fitting. The points represent the photoelectron signal intensities 
Is Of bicarbonate as a function of your bulk molar fraction X,,: The concentration variates from 
the different solutions pH. 
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CHAPTER 6 - CONCLUSION 


CO, and its presence in the atmosphere is one of the primary factor influencing 
global warming. Moreover, the presence of this molecule in aqueous solutions such as 
oceans is responsible for their acidification. On top of that, its presence in aerosols may 
increase the amount of acid rain. Since all these phenomena involve the adoption of CO, 
in aqueous solutions and this is a surface phenomenon. It seems to be of importance to 
study the presence of CO, and related deprotonated species in that region. Therefore, in 
this thesis, the XPS technique with microjet technique was employed to study this problem. 
The essential characteristics of the x-ray spectroscopy for the present study are the element 
specificity and the surface sensitivity. 


In this work, the XPS spectra used to comprehend the distribution of protonated 
molecules in the aqueous surface when carbon dioxide is dissolved in liquid water for 
different solution's acidity levels. The chemical shit between both deprotonated species was 
enough to identify different species. The spectrum at pH 7.0 was the only one which shows 
the presence of gas CO,. According to our experimental results, the surface enrichment 
model turned out to be excellent in explaining the observations as compared to the 
Hendersen-Hasselbalch model. With the help of this model, some crucial aspects of the 
surface distribution of the species study in this work as a function of pH were obtained. It 
was possible to conclude that both carbonate (HCO,”) and bicarbonate (CO?) are deprived 
of the surface. In particular, bicarbonate (CO?) is even more depleted from the surface than 
carbonate. The pK, value found for species at the surface is smaller than the bulk value. 
The explanation for this was discussed, and the temperature is one of the factors that will 
shift the pK,. The amount of change in the pK, is given by the Van Hoff equation. However, 
the expected value predicted by the equation is higher than the one at room temperature. 
Therefore one must rule out the temperature as a determinant factor. Another explanation 
for the observation is the presence of the third species, CO,, at low pH which steals intensity 
from the carbonate (HCO,”) which is not supported by the surface enrichment model. 


Another technique used in this work was NMR. Based on the chemical shift found in 
the NMR spectra, the molar fraction of bicarbonate and carbonate in bulk was obtained. The 
pK, in bulk was also determined. 


A comparison between the concentration variation of the two charged molecules, 
carbonate (HCO,) and bicarbonate (CO,?), as a function of pH for the surface and bulk 
regions, was carried out. A considerable change in the concentration of the two species 
distribution between the two areas was observed. For higher pH, the distinction is more 
accentuated. 


From our results, it was possible to conclude that both charged species are not prone 
to stay on the surface. The doubly charged species avoids the surfaces even more than the 
singly charged. 


Chapter 6 - Conclusion FER 


The Langmuir isotherm model is often used in the community to determine the 
surface concentration a species in aqueous solutions as a function of its mol. Still about the 
Langmuir isotherm model, in this work, it was proposed a new approach to obtain the Gibbs 
free energy. The variation of the concentration of the species of interest was monitored as 
a function of the pH. The surface concentration was obtained using the XPS signal intensity 
while the bulk was obtained experimentally by RMN data. This new proposed procedure 
allowed the evaluation of the Gibbs free energy for bicarbonate. A less satisfactory value 
was obtained for carbonate because this molecule participates in two chemical reactions, 
which is not taking into account by the model. 


The use of x-rays spectroscopy to study aqueous surfaces was a significant advance. 
Quantum mechanical simulations could help bring more information to the low pH region. 
Relative XPS intensity of the two species studied as compared to a third well knows species, 
such as formate, could give the needed knowledge of how much each species are present 
in the surface compared to a neutral species. The difference in surface propensity between 
a neutral species and singly change one may be much more significant than the observed 
differences between singly and doubly charged species. 


Chapter 6 - Conclusion 
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APPENDIX A ERROR ANALYSIS 


An important step after the data analysis is the error analysis, that will help determine 
the validity of our data and the Langmuir model. The area above the peaks are the interest 
of this work, which is proportional to the amount of a that certain molecule in the sample. 
Remembering that we have five samples at different acidity with three solutes in each 
sample. As explained in section 2.4, there is a change in concentration of the three species 
depending on the acidity of the medium. We can observe this variation in the measures by 
the intensity, see figure 5.1. 


The error analysis was perform from the intensity of each peak. The intensity is given 
by the number of counts over a period of time until the next energy step.This intensity is not 
absolute, there is a distribution of those counts. The fit gives us the most probable intensity 
from the parameters that we defined at the beginning of the analysis. Performing small 
variations in the intensity of each peak, we find maximum and minimum intensity values 
within the distribution and as consequence the maximum and minium areas that will be used 
to extracted the error bars. For this procedure the width at half height of the gaussian and 
the lorentzian are fixed as well as the energy position. 


— Background — Background 
Dispersal O Dispersal , 
— Fit curve - Maximum intensity — Fit curve - Minimum intensity 


C1s PE intensity [Arb. units] 
Cis PE tensity [Aro. units] 


T T T T T T 
292 294 296E 288 290 292 294 296E 
Binding energy (eV) Binding energy (eV) 
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Figure A.1: The (a)maximum and (b)minimum intensity of 9.8 pH. 
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288 290 


To better illustrate what was done we are going to take as example the bicarbonate 
from the 9.8 pH. Figures A.1 show the spectrum of that sample, in each figure we have the 
maximum and minimum intensity, respectively. We are interested in the normalized areas of 


each specie, that is, the normalized area of bicarbonate Area co; nom IS 
7: Area co; 
rea -morm 
HCO3 snorm Areaco, + Areancos + Areagoz- (A.0.1) 
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where Areaç,, is the area ofthe CO,, Area, co; is the bicarbonate area and Area o2- is 


HCO 


the carbonate area. In the case of figures A.1 we have the areas Area, coz max and Area, cos min 
for the maximum and minimum intensity of bicarbonate. The maximum and minimum 


normalized area are given by the combination 


Areâpco; 


mar 


Area = (A.0.2) 
< HCO, ;norm;mazx E . sa Are, N Pe 
: Areacos min + Area co; mar + Áreacoz- min 


Area pco; 


min (A.0.3) 


Area - n E 
= HCOs ;snorm;min o . 
+ Ar aco: 


Areacos mar + Area pico; 


«min «mar 


The same procedure is performed for the others species and for all samples. The 
error bars are given by the difference between those areas and the areas found by the 
fitting, that is 


Errormas = AreG CO; normmar - Are co; norm (4.0.4) 


Aro [a 
morm A CACO; morm,min (A.0.5) 


Errormin = Areâgco; 


Figure A.2 is a representation of the error bars calculated above. 


+ ENO yax 


+ ErrOr min 


Figure A.2: Representation of a data point with yours respective error bars. 
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